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THEORETICAL AND EXPERIMENTAL INVESTIGATION OF THE EFFECT OF TUNNEL WALLS 
ON THE FORCES ON AN OSCILLATING AIRFOIL IN TWO-DIMENSIONAL 

SUBSONIC COMPRESSIBLE FLOW 1 

By Harry L. Runyan, Donald S. Woolston, and A. Gerald Rainey 


SUMMARY 

This report presents a theoretical and experimental investi- 
gation of the effect of wind-tunnel walls on the air forces on an 
oscillating wing in two-dimensional subsonic compressible 
flow. A method of solving an integral eguation which relates 
the downwash on a wing to the unknown loading is given, and 
some comparisons are made between the theoretical results and 
the experimental results. A resonance condition , which was 
predicted by theory in a previous report (A 7 AC A Rep. 1150), 
is shown experimentally to exist. In addition, application of 
the analysis is made to a number of examples in order to illus- 
trate the influence of walls due to variations in frequency of 
oscillation, Mach number, and ratio of tunnel height to wing 
semichord. 

INTRODUCTION 

In the evaluation of results obtained by measurement of 
the forces on a wing in a wind tunnel, the question of the 
effect of the tunnel walls arises. In the case of steady flow 
the problem has been extensively investigated and, in gen- 
eral, relatively simple factors have been determined which 
can be used to modify measurements of the forces on a wing 
in a tunnel to correspond to free-air conditions. However, 
the corresponding problem of the effect of walls on an os- 
cillating airfoil has received relatively little attention, par- 
ticularly in the case of compressible flow. The present 
report concerns the wall effects in the oscillating case and 
treats the problem in two-dimensional subsonic compressible 
flow. 

In incompressible flow, theoretical treatments of wall 
effects on oscillating wings have been made by several in- 
vestigators and reported in references 1, 2, and 3. These 
investigators have shown generally that the tunnel-wall 
effects are a maximum for some small values of the reduced 
frequency and that the wall effects become negligible as the 
reduced frequency is increased. Extension of the theoretical 
treatment of the problem to include the effects of compressi- 
bility of the fluid has been reported in reference 4. In this 


reference, it is shown that, in addition to the large effect 
noted at low values of the reduced frequency, under certain 
conditions, large effects of the walls may be encountered at 
higher values of the reduced frequency. These effects are 
due to an acoustic resonant phenomenon which occurs when 
a disturbance from the oscillating wing is reflected from the 
tunnel wall back to the wing with such a phase relationship 
that it reinforces a succeeding disturbance. 

In reference 4, the problem was expressed as an integral 
equation which relates the known downwash distribution 
over the airfoil to the unknown lift distribution. One pur- 
pose of the present report is to discuss further the integral 
equation and to demonstrate a method of solving it. A 
second purpose is to present some results showing wall effects 
calculated by this procedure and, in some cases, to compare 
the calculated results with experimental results. This phase 
of the investigation is given in three parts: (1) A comparison 
between analytically and experimentally determined values 
for the lift and moment on a wing oscillating in pitch at 
several subsonic Mach numbers; (2) an analytical study of 
the effects of a variation in Mach number for a constant ratio 
of tunnel height to wing semichord; and (3) an analytical 
study of the effects of a variation in the ratio of tunnel height 
to wing semichord. Portions of this material have been 
reported previously in reference 5 and are included in the 
present report in order to provide a more extensive and 
unified presentation. 

As a check, the integral equation for the downwash on a 
wing oscillating between walls in a compressible medium is 
reduced to the zero-frequency condition and is given in the 
appendix. The resulting expression is in agreement with 
steady-state results. 

The calculation procedure and the results contained in this 
report are of significance for such problems as the experi- 
mental measurement of the forces on an oscillating airfoil, 
the determination of wing-flutter characteristics in wind 
tunnels, and also in certain possible types of flutter of airfoils 
in cascade. 


1 Supersedes NACA Technical Note 3416 by Harry L. Runyan, Donald S. Woolston, and A. Gerald Rainey, 1955. 
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SYMBOLS 


a 


A 


n 


b 

h 

71 

II 

k 

L(x 0 ),L(do) 

La 


M a 

Mn 

M 

Mk II 
P 2tt/3 


velocity of sound, ft/sec 
coefficients in series expression for 
lift distribution (eq. (16)), where 
n= 0, 1,2,... 
wing semi chord, ft 
displacement of wing in vertical 
translation, ft 
height of tunnel, ft 
height of tunnel referred lo wing 
semi chord 

Hankcl functions of the second kind 
reduced-frequency parameter, bco/U 
kernel of integral equation 
lift distribution, lb/ft /unit span 
aerodynamic lift force per unit span 
due to pitch 

aerodynamic lift force per unit span 
due to translation 

aerodynamic moment per unit span 
due to pitch 

aerodynamic moment per unit span 
due to translation 
Mach number, Ufa 


R n ^^{x~x Q ) 2 A^\nH) 2 , where ^=1,2,3, . . . 


U 

w{x) 

X a 


x t xo,y£ 
z=Jc(x—x 0 ) 


stream velocity in chordwise direc- 
tion, ft /sec 

vertical induced velocity (perpen- 
dicular to chord), ft/sec 
axis of rotation measured from mid- 
chord, positive rearward, based on 
semi chord 

Cartesian coordinates 


a 


angular displacement of wing in 
pitch, radians 


18 - Vl -M 2 


€ = 


k— Ail 

ffl 


6=— cos k 
d 0 ~— cos _ k 0 


p 

K 

<t>L u 

hi a 

4>M h 

CO 

w res 


fluid density, slugs/cu ft 
phase angle between lift force and 
position of pitching wing, deg 
phase angle between lift force and 
position of translating wing, deg 
phase angle between moment and 
position of pitching wing, deg 
phase angle between moment and 
position of translating wing, deg 
circular frequency of oscillation, 
radians/sec 

circular frequency at resonance, 
radians/sec 


A p pressure difference between upper 

and lower surface, lb/sq ft 
Primed quantities refer to a wing in free air. 

ANALYTICAL INVESTIGATION 

This section is concerned with the development of a method 
for solving the integral equation, originally derived in refer- 
ence 4, which relates the downwash to the loading on an 
oscillating wing. The basic integral equation and its kernel 
is given by equations (1) and (2). Reduction of the kernel 
is made in equations (3) to (10). Alternative series expres- 
sions for the kernel which are suitable for numerical computa- 
tion are given by equations (11) to (15). The loading on the 
wing is given by equation (16), the downwash expression by 
equations (18) and (19), and finally the lift and moment 
expressions by equations (20). 

THE INTEGRAL EQUATION AND ITS KERNEL FUNCTION 

The integral equation. The integral equation of reference 
4 for the vertical velocity or downwash of an oscillating airfoil 
between plane walls may be written as 

w(x) = J L( Xo ) [K(M, z) + K (M, z,II)]dx 0 (1 ) 


where w(x) is the known vertical velocity (or known motion 
of the wing) and L(x 0 ) is the unknown lift distribution or the 
local strength of a distribution of oscillat ing pressure doublets. 
The functions within the brackets comprise the kernel func- 
tion of the integral equation and appear formally as 


K(M,z)=jL r Jim e ~ mz ~ Xo) 
4 p/C 0 


K(M,z,H)=~ lim " 2 2 (-1)* 

4p/C y— >0 J — oo v = I 


(2a) 


J« S 2 


dif 


Ijt II o <2 > VC+0 2 (y-«^) 2 ] da (2b) 


The first function K(M , z) corresponds to the kernel for the 
free-air condition as given by Possio (ref. 6). The second 
function K(M, z, II), containing the infinite summation, is 
the additional part of the kernel arising from the effect 
of the walls. Physically, a kernel function represents the 
contribution to the vertical velocity at a field point due to a 
pulsating pressure doublet of unit strength located at any 
other point in the field. For the particular case represented 
by equations (2), the kernel function gives the vertical veloc- 
ity in the plane of a wing located in the center of the tunnel. 
The expression K(M, z) gives the downwash of a doublet in 
die plane of the wing, whereas the expression K(M, z, H) 
gives the downwash due to the system of images which mathe- 
matically represents the walls. 

Reduction of the kernel function — The integrals contained 
in the expressions for the kernel function inequations (2) are 
improper because they have an infinite limit and also be- 
cause. at certain points, the integrands become singular. 
This section is concerned with the reduction of these integrals 
to a form more amenable to computation. 
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By making use of the fact that the Hankel functions in equations (2) satisfy the identity 

U // o (2 > (~ VFw)-^ ^o <2> VFW) 

there is obtained for the downwash 

«<*>-$ JS «?. J ‘ w-' ('>’/_? ^'s? »•* (f V!W) rf£+ 

■ ' #o (2) #W) #+ 


(3) 


m 2 p r*-*» 

p- 


2±(~1 )* {/3 2 J‘ [y\r d 2 ('/-«//) 2 ]r/H- 

^ J_ r “ r "e^ W> Vl 2 +)3 2 (2/-nH) 2 ] d?}) 


M 2 ¥ 


(4) 


The integrals of equation (4) that contain partial derivatives of Hankel functions can be integrated twice by parts to 
obtain 


f** ^ Ho <2> V€*+/S*(lf-nfl)*] « - • 




P 2 




ik s^o) 


V(*-*o) 2 +/5 2 (2f-«fl) 2 ]-^ J_Y° e ? #o (2) [^T <ft 


(5) 


The last integral of equation (5) may be written in two parts 
as 

ik 

Zj (-1)” f^e* 2 W> [— Vr+Xlf-nfl)*] 


j: 




8 - 


S(-lW 2 >[^VFh8% 


</?+ 


-«H) 2 ] 


£ s (-l)”-^o (2) -Je+Fiy-nH) 2 ] di 


(6) 


The first integral on the right-hand side of equation (6) will 
be left temporarily in integral form and will be treated in the 
following section. (See evaluation of £ 3 following eq. (13).) 

The second integral on the right-hand side of equation (6) 
has not been integrated in closed form; however, in wind- 
tunnel problems it can be handled conveniently by approxi- 
mate methods. (An alternative means of treating this inte- 
gral, which avoids the approximation but is somewhat more 
tedious, will be indicated in the discussion following eq. 
14(c).) A practical assumption which is often made in the 
analysis of the effect of wind-tunnel walls is that the tunnel 
height is considered large compared with the wing semichord. 
With this assumption the argument of the Hankel function in 
equation (6) can be written as (in the limit as ?/— > 0) 


may be replaced by the vertical distance finll of the image 
above the airfoil. Of course, this approximation does not 
hold for Mach numbers close to or equal to unity. With this 
approximation, the second integral of equation (6) can be 
expressed as 

J Cx-x o fe oo FMJr , 1 

e* S(-l)'7/„< 2 > 1V p^e+f3 2 (y-nI/r 
'0 11 = 1 L P J 

— rn 


= £ (-l)'7/o (2 > ( M ^ // )J 1 X 


(7) 


and these equations may be used to express equation (4) as 
cob f 1 

w(x)= ^ , J ^ L(x 0 ) [ K(M, z) + K{M , z,H)\ dx 0 (8) 


where 


K(M, z)=^ e-*‘ | 


+ 

(9) 


and 


^aV+/3 2 (nny=^ Pull 


vG is) 


, . Mk 
+1 ~-0- nH 


provided that ■ 1 . 


This approximation implies that the airfoil images, and, 
particularly the closest image n=l, are a sufficient distance 
from the airfoil so that the actual distance fi 2 (riH) 2 


K (M ,z,H)= e ~'p [-e^ ± ( — 1)"//„ <2) (nli n ) + 

P 2 (e*-l)± i-iyi-u 2 ' (~^)+ 

ik £ e~¥ II 0 m [^£W) ! ] di + 

J £ (-i)» ? tf M 7/, (2 >wl (10) 

71 = 1 J 
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in which use lias been made of 


kM 



Mz 


v= 




li 0 =\x— X 0 | R n = y ! (x — Zo) 2 + 0 \nH ) 2 


Equation (8), together with the definition of equations (9) 
and (10), permits the determination of the effect of tunnel 
walls on a lift distribution L( 6 0 ) for a given downwash distri- 
bution w(x). The integral equation for the case of no tunnel 
walls checks the results of Possio (ref. 6). For the case with 
walls and for the limiting steady -flow case of zero frequency, 
it is possible to obtain a mathematical check with some 
existing results; this check is shown in the appendix. 

Alternative series expressions for kernel. Although the 
form of the kernel given by equation (10), could 

be used for calculation, alternative series which are more 
highly convergent may be used and are given in this section. 

The kernel K(M,z,H) is the sum of four infinite series 
which can be written as 


k(m,z,h)=- 2 ~ (aSi+CiSt+asz+CtSd ai) 

where the S n ’ s denote the indicated infinite summations of 
equation (10) and the C n ’ s the respective multipliers. 

Series Si and S 2 of equation (11) may be put in a more 
rapidly convergent form according to Infeld, Smith, and 
Chien (ref. 7). When the variables <p and e are introduced, 
where 

MkH 
V= 


and 


2t rfi 


\x—x^ 

m 


the series Si and S 2 can be written as 

s^±(-\r iio (2) (M 

77 = 1 

= (- D" W 2 > (2x V V« 2 +n 2 ) 


1 r 2 ie ~ 7re '^ 1— 4 
~ 2w L \ 1 — br 


4 p- co 

+ 2 iS 


n = I V(2»+l) 2 — 4^ 2 


/(?»- l) J -4p2 


' — it Ho <2) (27r ep) 

V(2ra-l) 2 -4p 2 


] 


( 12 ) 


and 

s 2 =s(-i)"//o <2 > CT) 

= S(- 1)* // o <2) (27 rnp) 

77 = 1 

where Euler’s constant 7=0.577215. 


(13) 


Series S 3 may be evaluated by utilizing the expression for 
Si (eq. (12)) and integrating the resulting expression to 
obtain 


S 3 =r e 'k ± (-1)" 7/o <2) (14a) 


*--y. 

2ifS 


» oo ik 

( - ii 




IVI 77=0 




di 


(14b) 


v 0 , 1 +/S , 

S 3 = ~,k l0Sc M + 


210 


2 


Al n-o ( 7T \~ 


M 

k \ 2 P 2 H 


(m) [(2n+1)2 - 4 ^+(|) _ 

k/P 2 


TVp 

M 


V(2^+l) 2 — ^V 2 


(14c) 


i t is of interest to note that series S 3 may be employed in 
an alternative means of integrating equation (7). For 
application to wind-tunnel problems, where the ratio of 
tunnel height to wing semichord is small, or in application 
to cascade problems, the approximation employed in integrat- 
ing equation (7) becomes less valid. It is possible to avoid 
the use of the approximation by writing the integral of 
equation (7) in a form which is identical to that of equations 
(14a) and (14b) with the exception of the upper limit. The 
integral containing the Hankel function can be evaluated 
by employing the tables of Schwarz (ref. 8). The second 
integral, containing only an exponential term, can be in- 
tegrated in closed form, as was done to obtain equation (14c). 

Series Si may be evaluated in a direct manner by employ- 
ing tables of the Hankel function and by using for large 
values of the argument the approximation 

2 (15) 
V TT/dUn 

With the aid of series S u So, S 3 , and S i} the kernel K(M,z,II ) 
may be evaluated. 


METHOD OF SOLUTION 

A method of using equation (8) to determine the aero- 
dynamic forces on a wing oscillating in the presence of plane 
walls is now discussed. The method under consideration is 
one of collocation similar to that used by Possio (ref. 6) 
and Frazer (ref. 9) for the case of no walls. The approach 
involves the assumption of an appropriate series expression 
for the lift distribution, substitution of this series in the 
integral equation for the downwash, and calculation of the 
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downwash at arbitrarily selected points on the chord (con- 
trol points). Thus equation (8) is reduced to a set of simul- 
taneous equations, the unknowns of which are the coefficients 
of the assumed expression for the loading. 

Expression for the loading. -The expression which is 
assumed for the lift distribution is a trigonometric series 
expansion which satisfies the Kutta condition at the trailing 
edge and which has the proper type of singularity at the 
leading edge. This expression is 

— YtY==Ao cot A n sin 7160= L(6o) (16) 

pU 2 n = i 

where ar 0 = — cos 0 O and the A n ’ s are unknown coefficients 
to be determined in accordance with the downwash w{x) } 
which is known from the motion of the wing. It is desir- 
able to rewrite equation (8) in terms of the variable 0 O 
as follows: 

?^(x)=C/A* L(do)K(M,z)s\n6 o d0o+ j i(0 o )/f(M,2,fZ)sin0 o d0oJ 

(17) 


The first integral on the right-hand side of equation (17) 
is the integral expression first derived by Possio (ref. 6) for 
the condition of no walls. Its solution has been treated by 
several investigators (see, for example, ref. 9) and will not 
be discussed herein. It can be expressed entirely in terms 
of the unknown coefficients A n of equation (16). The sec- 
ond integral of equation (17) may be evaluated by the use 
of equations (12), (13), (14), and (15). 

Determination of the aerodynamic forces. The integrals 
of equation (17) are determined for a selected number of 
control points and equated to the expression for the down- 
wash. The expression relating the downwash to the motion 
of a wing translating ( h ) and pitching (a) about an axis 
located at x a is 

w(x) =h-\-Uot-\- b(x — Xa)oc (18) 

or, with the assumption of harmonic motion, 

^~-=ik -|+[1 + ik{x—x a )\<x (19) 


Equation (19) is used to calculate w(x) for values of x ap- 
propriate to each of the selected control points. A set of 
simultaneous equations can then be written, the number of 
which corresponds to the number of control points employed 
and (conveniently) to the number of terms retained in the 
series for L(0 O ). The unknown coefficients may now be 
determined by solving these simultaneous equations. The 
total lift and moment about the midchord are given in terms 
of the coefficients A n through the relations 


-L a= \ 

TrpbU 2 2 



> 


M a _1 
t rpb 2 U 2 8 



( 20 ) 


Effect of the number of control points considered. — An 
investigation was made of the number of terms of the series for 
the lift distribution (eq. (16)) and thus of the number of 
control points required to obtain satisfactory accuracy. 
Calculations were performed for a particular case by in- 
creasing the number of control points and the number of 
terms of the loading series until the solutions were in reason- 
able agreement. For the case considered, three terms of 
the series for the lift and three control points at the quarter-, 
half-, and three-quarter-chord positions gave satisfactory 
results. The consideration of two additional control points 
at the leading and trailing edges, together with two addi- 
tional terms of the lift series, made no significant change in 
the results. For high values of the reduced-frequency 
parameter k } the use of additional control points might be 
necessary. 

The procedure just discussed involves consideration of a 
continuous distribution of pressure doublets over the chord. 
Calculations requiring much less computing can be made by 
considering the chord wise loading to be concentrated in a 
single doublet located at the quarter chord and by satisfying 
the downwash at the three-quarter chord. In the case of 
the lift, this approach has been found to give fairly good 
agreement with the results of the more elaborate calcula- 
tions except in the vicinity of the resonant frequency. 

THE ANALYTICALLY INDICATED RESONANCE PHENOMENON 

Two-dimensional tunnel. — By examination of equations 
(12) and (13), it may be seen that the series Si and S 2 become 
infinite when 

4p 2 =(2n-l) 2 

or where 

— =T0(2n-l) (71=1,2,3, . . .) (21) 

a 

At these critical values of the frequency parameter, the 
expression for the kernel K(M,z,H) (eq. (11)) becomes 
infinite for all values of x. Physically, this condition 
represents a resonance in the tunnel involving a transverse 
oscillation of the moving air between the walls. 

The fundamental or smallest critical values of coH/a 
corresponding to n= 1 in equation (21) are shown plotted as 
functions of Mach number M in figure 1. Equation (21) 
and figure 1 show that finite values of the critical frequency 
exist for the condition M= 0, 17=0, and . These 

conditions correspond to a compressible fluid at zero velocity 
in the tunnel. As the Mach number is increased, the 
critical-frequency parameter decreases rapidly and becomes 
zero at a Mach number of unity. 

As indicated by equation (1), the product of the lift and 
the kernel function must remain equal to the vertical 
velocity over the wing; this velocity is defined by the motion 
of the wing and remains finite. The product of the lift and 
the kernel function can remain finite only if the lift ap- 
proaches zero as the kernel becomes infinite. This condition 
in the tunnel is analogous to t he well-known case of a simple 
undamped-spring-mass system for which, at the resonant 
frequency, theory predicts an infinite deflection of the mass 
occurring even with a forcing function of small amplitude. 


Resonant frequency parameter, 
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Figure 1. — Fundamental values of resonant frequency parameter — as a function of Mach number. 
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Circular tunnel. -A resonance can also be demonstrated for 
the infinite circular tunnel. The nature of the boundary- 
value problem, for this case, makes it possible to separate 
variables; therefore, the governing partial-differential equa- 
tion can be reduced to Bessel’s equation. (See, for instance, 
ref. 10.) The resonant frequencies are then found as the 
roots of the equation 



where J n represents the Bessel function of the first kind, D 
is the tunnel diameter, and p n is the root of the equation 

Jn(Pn) = 0 

Values for p n for the first several modes are p n = 1.84, 3.05, 
and 4.17. Note that, for a circular tunnel having a diameter 
equal to the height of a plane tunnel, the fundamental 
frequency is 3.68/7t= 1.17 higher than resonant frequency 
in the plane tunnel discussed in this report. 

EXPERIMENTAL INVESTIGATION 

WIND TUNNEL 

The experimental part of the investigation of the effect 
of tunnel walls on the forces acting on an oscillating airfoil 
was conducted in the Langley 2- by 4-foot flutter research 
tunnel. For these tests, a rectangular test section having 
dimensions of 2 feet by 3.8 feet was used. This tunnel is of 
the closed-throat, single-return type and employs either 
air or Freon-12 as a testing medium at pressures from 1 
atmosphere down to about % atmosphere. 

It has been shown previously that the resonant frequency 
varies directly as the speed of sound. Inasmuch as Freon- 12 
has a speed of sound equal to about one-half that of air, the 
experiments to be discussed were conducted in Freon-12 
so that the resonant frequency could be surveyed within 
he frequency limitations of the equipment. 

MODEL AND OSCILLATING MECHANISM 

Figure 2 is a schematic drawing of the test section with 
the model and oscillating mechanism installed. The model 
had a chord of 1 foot and an NACA 65-010 airfoil section; 
it completely spanned the 2-foot dimension of the test 
section. The gaps between the model and the tunnel wall 
were sealed by end plates which rotated with the model. 
The model, driven symmetrically from both ends, was 
oscillated in pitch about the midchord by a direct-drive 
eccentric-cam system powered by an induction motor with 
variable frequency supply. 

INSTRUMENTATION 

The lift and moment on the wing were obtained by 
electrical integration of the outputs of 12 model 49-TP 



Bearing J L" 875 8 4 

Figure 2. — Schematic drawing of test section with model and oscil- 
lating mechanism installed. 

NACA miniature electrical pressure gages. The pressure 
gages, which are described in considerable detail in reference 
11, were located at the center of the span at 2.5, 5, 10, 15, 
20, 30, 40, 50, 60, 70, 80, and 90 percent of the chord. Each 
gage was arranged to indicate the difference in pressure 
between orifices on the upper and lower surfaces. Electrical 
integration techniques used in these experiments are dis- 
cussed in reference 12. The so-called square-wave method 
of weighting was used; that is, the pressure indicated by 
each gage was assumed to represent the pressure acting 
over a portion of the chord extending one-half the distance 
to the next gage both forward and rearward. For example, 
the fraction of the chord assigned to the first gage was 3.75 
percent and to the sixth gage was 10 percent. Some of the 
implications of this method of integration will he discussed 
in a subsequent section. 

The angular displacement at the midspan position was 
indicated by resistance- wire strain gages attached to a 
torque rod running through the center of the hollow wing. 
One end of the torque rod was fixed to the center of the 
wing and the other end was fixed to the tunnel wall. 

A schematic diagram of the instrumentation is shown in 
figure 3. The magnitude of the vector representing the fun- 
damental component of lift or moment and angular displace- 
ment was indicated on an alternating-current vacuum-tube 
voltmeter attached to the output of a variable-frequency, 
narrow-pass-band filter. In essence', the filter performed 
the function of a Fourier analysis in that both random com- 
ponents and higher harmonics were removed from the signal. 
In order to measure the phase angle between lift or moment 
and the angular displacement, the output of the filter was 
fed into a pulse-shaping circuit designed to convert the sinus- 
oidal signals into pulses corresponding in time to the “cross- 
over’* points of the original signal. The pulses were then fed 
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Fkii re 3. — Schematic diagram of instrumentation. 


into an electronic chronograph that accurately indicated the 
time interval between the leading pulse which started the 
chronograph and the lagging pulse which stopped it. The 
ratio of this time interval to the period of oscillation, when 
multiplied by 360°, yields the phase angle in degrees. The 
period and frequency of oscillation were determined by start- 
ing and stopping the chronograph with the angular-dis- 
placement signal. In order to minimize the effects of small 
differences in components between the two circuits, a “tare” 
switch was provided which fed a single signal (the angular 
displacement) through both circuits. The resulting phase 
angle represented the phase shift introduced by the filters 
and pulse-shaping circuits. 

TEST CONDITIONS 

The Mach number of the tests was varied from A/= 0.35 to 
]l/=0.7 and the Reynolds number was held constant at about 
5X10 6 by varying the density. The frequency of oscilla- 
tion was varied from 0 to 60 cycles per second, and the mag- 
nitude of angular displacement was about 1.2° except for 
some lift data at M= 0.71 which was obtained at an angular 
displacement of about 2.4°. 

DISCUSSION OF RESULTS 

The theory and calculation procedure and the experimental 
technique discussed previously for the determination of the 
forces acting on a wing oscillating between walls have been 


applied to a number of specific examples. The investigation 
has been divided into three parts: (1) A comparison is 
made of analytical and experimental results obtained for 
the lift and moment on a pitching wing for several subsonic 
Mach numbers, (2) theoretical results for the effects of a 
variation in Mach number at constant tunnel height are 
given for a pitching wing and also for a wing undergoing 
vertical translation, and (3) theoretical results for the effects 
of a variation in the ratio of tunnel height to wing semichord 
are presented for particular values of Mach number. 

COMPARISON OF THEORY ANI) EXPERIMENT 

In figure 4 a comparison is made of analytical and experi- 
mental results for a wing oscillating in pitch about its mid- 
chord. Figures 4 (a), 4 (1)), 4 (c), and 4 (d) apply, respec- 
tively, to Mach numbers of 0.35, 0.5, 0.6, and 0.7. The 
results apply to a ratio of tunnel height to wing semichord 
II of 7.60. 

The plots on the left-hand side of each figure show the 
magnitudes of the forces and moments as a function of the 
frequency of oscillation, whereas those on the right-hand 
side show the corresponding phase angles. The magnitudes 
are presented as ordinates in the form of ratios \L a /L a '\ and 
\MJM a '\. In these ratios, the quantities L a and M a are, 
respectively, the lift force and the moment on a wing in a 
tunnel; LJ and MJ are the theoretical lift and the theoreti- 
cal moment on a wing in free air. The effect of the tunnel 
walls appears, therefore, as a deviation from unity of the 
ratios \L a /L a '\ and \MJM a f \ when L a and M a are the theo- 
retically derived forces and moments. When L a and M a 
represent the experimental forces and moments, the devia- 
tion from unity may not be completely attributed to the 
effect of tunnel walls because such factors as airfoil thick- 
ness and viscosity may cause deviation from the elementary 
theory. The abscissa in the figures is the ratio of the fre- 
quency of the pitching oscillation to a frequency calculated 
for the resonant condition. 

Excellent agreement between theory and experiment is 
obtained for the phase angles, in most cases, for both the 
lift and the moment. Quantitatively, however, the agree- 
ment between theory and experiment for the magnitudes 
of the forces is not as good, although very similar trends 
are demonstrated; in most cases, a systematic difference 
appears. Some possible sources of the differences between 
theory and experiment are discussed in the following section. 

Examination of figure 4 reveals that the theory predicted 
the resonant frequency very well. In all cases, the mini- 
mum lift and moment were found to lie very close to the 
analytically indicated resonant frequency. Theoretically, 
the lift and moment reduce to zero at the resonant condition. 
Under actual conditions, such as finite tunnel length, trans- 
mission of energy through the walls, nonlinearities at higher 
amplitudes, and turbulence in the flow that gives rise to 
damping, pure resonance is unobtainable. However, it 
may be seen by examining figure 4 (d) that the lift and 
moment were reduced to 20 percent of the values away 
from resonance. 
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Pigi're 4. — Comparison of theoretical and experimental results for the magnitudes and phase angles of the lift and moment of a pitching wing. 

Height-semichord ratio H = 7. (SO. 
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(b) M = 0.5. 

Fk; i rk 4. — Continued. 
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(c) A/ =0.6. 
Figure 4. — -Continued. 
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(d) M =0.7. 
Figure 4. — Concluded. 
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REMARKS ON SOME DIFFERENCES BETWEEN THEORY AND EXPERIMENT 

In this section, some limitations and possible reasons for 
differences between theory and experiment are listed and 
discussed. 

In the comparison between theory and experiment shown 
in figure 4, an almost constant difference of 10 to 15 percent 
between the magnitudes is to be noted, whereas the phase 
angles are in good agreement. Of the several possible 
reasons for these differences between theory and experiment, 
perhaps the more important ones are airfoil thickness, 
Reynolds number, finite tunnel length, transmission of 
energy through walls, dissipation of the pressure waves due 
to turbulence, and integration procedures. The effect of 
all these possibilities is not known for the oscillatory case. 
In the steady-state case, however, it is known that the 
effect of increasing thickness is to increase the slope of the 
lift curve. 

In the considerations of the analytical integration, a 
collocation scheme was used to solve the integral equation. 
In general, the three collocation points used were found to be 
satisfactory as pointed out previously. 

Twelve pressure cells were used for the experimental 
integration of the forces. A stepwise integration pro- 
cedure was employed; that is, the pressure as record- 
ed by the pressure cell was multiplied by an area of 
the wing over which it is assumed that the cell will give 
an average pressure. This procedure gives good results 
except possibly at the leading edge where the pressure varia- 
tion is very great. Theoretically the pressure approaches 

infinity at the leading edge ^as 1 ^ and experimentally it is 

found to be very large. As a matter of fact, if a theoretical 

distribution of pressure is assumed to be cot 1 ' as 

2 V 1 +x 

shown in the following sketch and this curve is integrated 



in the same manner as the experimental curve was inte- 
grated (that is, by calculating the ordinate at the same 
value of the abscissa at which the pressure cells were located 


for the experiment), it is found that the area as determined 
by the approximate method is 8 percent less than the area 

as determined by integrating cot 0 in closed form. It is 

apparent that the neglected area (shaded) can he appreciable. 
In the actual experiment, in which a fairly thick airfoil was 
used, the neglected area would probably be smaller but 
could perhaps contribute to the almost constant difference 
in magnitudes of lift and moment between the theoretical 
and experimental results. The effect of this integration 
difference on the phase angles, which have been shown 
to be in good agreement, would not be as pronounced. 

EFFECT OF A VARIATION IN MACH NUMBER AT CONSTANT TUNNEL 

HEIGHT 

An analytical investigation has been made of the effects of a 
variation in Mach number at constant tunnel height on the 
forces on an oscillating wing. Some of the results of the 
previous section are employed together with results of 
additional calculations. The magnitudes and associated 
phase angles of both the lift force and the moment have been 
determined for a pitching wing and also for a wing under- 
going vertical translation. Calculations have been made 
for a constant value of the height-semichord ratio II of 
7.60 and for Mach numbers of 0.3, 0.5, 0.7, and 0.8. Re- 
sults of the calculations are shown in figures 5 and 6 for the 
lift and in figures 7 and 8 for the moment. 

The magnitudes of the forces and moments are presented, 
as in the previous section, in the form of ratios: |L a /Z a 'j, 
|M«/M a '|, \L h /L h '\j and \M hl /M h '\. The phase angles re- 
lated to these ratios are presented as a difference oeiwueii a 
wing in a tunnel and a wing in free air. The magnitudes 
and phase angles are plotted against a frequency parameter 
c oH/a, where w is the circular frequency of oscillation of the 

wing, H is the height of the tunnel, and a is the velocity of 
sound. At a particular value of the frequency parameter, a 
progressively larger effect of the walls is indicated as the 
Mach number increases. At all Mach numbers, the lift 
is reduced to zero at the resonant frequency. The dip in 
the curves against frequency ratio, which appears to be 
characteristic of the low Mach number cases, gradually 
disappears as the Mach number is increased. As in the 
case of the magnitude of the lift, the effect of the walls on 
the phase angle increases as the resonant frequency of the 
tunnel is approached and also as the Mach number is 
increased. 

The magnitude of the moment about the midchord is 
shown in figures 7 and 8. The curves have the same shape 
as the lift-ratio curves and again decrease to zero at the 
resonant condition. The corresponding phase angles are 
shown in these figures. Note that in figures 5 to 8 only 
slight differences appear between results for pitch and those 
for translation. 


Phase angle, <£ -<A ‘,deg 
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Figure 5. Effect of variation in Mach number at constant height-semichord ratio 11 = 7.(>0 on the magnitude and phase angle of the lift of a 

pitching wing. 


Lift ratio, 
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Figure 6. — Effect of variation in Mach number at constant height-semichord ratio H= 7.00 on the magnitude and phase angle of the lift of a 

translating wing. 
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Figure 7. — Effect of variation in Mach number at constant height-semichord ratio //= 7.00 on the magnitude and phase angle of the moment of 

pitching wing. 
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Figure 8.— Effect of variation in Mach number at constant height-semichord ratio 11 = 7.60 on the magnitude and phase angle of the moment of a 

translating wing. 
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EFFECT OF VARIATION IN RATIO OF TUNNEL HEIGHT TO WING SEMICHORD 

The effect on the lift-force ratio and on the associated 
phase angles of varying the ratio of tunnel height to wing 
semichord is illustrated in figure 9. The results presented 
in figures 7 and 8 have been based on the consideration of a 
distribution of pressure doublets over the chord of the airfoil 
and on satisfying the downwash condition at three chordwise 
stations. Results for figure 9 have been obtained by the 
simplified procedure of concentrating the loading in a single 
doublet at the quarter chord and satisfying the downwash 
condition at the three-quarter chord. This approach gives 
fairly good agreement with the results of the more elaborate 
procedure except near the resonant frequency. 

Calculations have again been made for an airfoil oscillating 
in pitch about its midchord for values of the height-semi- 
ehord ratio 77 of 7.60, 16, and 32 at M= 0.3 and 0.8. In 
figure 9 (a) the lift ratio \L a fL a f \ is plotted against the re- 
duced-frequency parameter k = bco/U. Plots for both Mach 
numbers are made to the same scale for ease of comparison. 
It is again apparent that the effect of reducing the Mach 
number is to reduce the effect of the tunnel walls and to 
raise the value of the critical frequency at which resonance 
can occur for a given tunnel. For example, for 77=7.60 , 
at M= 0.8, the critical value of k is 0.30, whereas for M= 0.3, 


the critical value of k is increased to 1.31. Also, as was to 
be expected, increasing the height of the tunnel has a marked 
effect in reducing the influence of the tunnel walls for most 
of the frequency range. However, the critical frequency 
is reduced for increased tunnel height so that in the large 
tunnel the range of k below the fundamental resonance 
becomes smaller. This reduction in frequency would seem 
to be a disadvantage of the larger tunnels. However, the 
second branch of the curve for 77=32 at M= 0.3 shows that, 
for frequencies between the first and second resonant points, 
the effect of the walls on the magnitude of the lift is no 
greater than for frequencies below first resonance. Note 
also that the approach to resonance is quite abrupt; con- 
sequently, only a small range of frequencies very close to 
resonance would be critical and tests could then be conducted 
between the critical frequencies. 

In figure 9 (b) the phase angles, associated with the results 
of figure 9 (a), for both the wing in a tunnel and for a wing 
in free air are shown as a function of reduced-frequency 
parameter k for values of height-semichord ratio 77 of 7.60, 
16, and 32 at M= 0.3 and 0.8. 

At M= 0.3, the effect of walls on the phase angle is gener- 
ally very small; at M= 0.8, the effect is small at low fre- 
quencies but increases greatly as the critical value of k is 



(a) Magnitude of lift. 

Figure 9. — Effect of a variation in tunnel height on the lift of a pitching wing for M=0.3 and M=0.8. 
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(b) Phase angle of the lift. 
Figure !). — Concluded. 


approached. As the resonant frequency is approached, 
results at both Mach numbers show that the phase angle 
increases in negative value and appears to approach— 90°. 
As the resonant condition is exceeded, there is a sudden shift 
in phase angle. This change is similar to that found for the 
oscillation of a simple undamped spring-mass system where 
an abrupt change in phase angle of 180° is found as the 
resonant frequency is exceeded. Because of the complexity 
of the kernel, which involves an infinite series of Hankel 
functions, the phase angle at resonance has not been 
determined. 

In figure 10, the plot from figure 9 (a) for M= 0.3 and 
11= 7.60 is compared with some results of Reissner (ref. 2) 
for the effect of walls on the lift-force ratio in incompressible 
flow. This curve for M= 0 significantly duplicates the rather 
large wall effect at low values of k which is noted for the 
M= 0.3 result. (At low values of k y the curves for the two 
different Mach numbers are almost coincident probably be- 
cause of the fact that the slightly lower height semichord 
ratio H= 6.28 used by Reissner counteracts the effect of de- 
crease in Mach number.) For values of k greater than 0.5, 



Figure 10. — Comparison of lift ratio for M = 0 and M= 0.3. 
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the curves for the two Mach numbers separate; the lift-force 
ratio at M= 0 approaches unity and exhibits no effects of 
resonance because the resonant phenomenon arises only from 
the effects of compressibility. 

CONCLUDING REMARKS 

This report has dealt with the problem of an airfoil oscil- 
lating between plane walls in subsonic compressible flow. It 
constitutes a continuation of the work initiated in NACA 
Report 1 150 in that a method of solving the integral equation 
is presented and some experimental results are compared 
with theory. 

The comparison between theory and experiment for the 
phase angles between lift force or moment and position is 
shown to be very good, whereas the comparison between 
theory and experiment for the magnitudes of the lift and 
moment is not as good; however, the trends are all accurately 
predicted. In all cases the resonant frequency was accu- 
rately predicted. The cause of the apparent discrepancy be- 
tween the theoretical and experimental lift and moment may 
be attributed to several factors such as dissipation of the 
pressure waves due to turbulence of the air flow and trans- 
mission of the energy through the tunnel walls. In addition, 


the theoretical work was based on the concept of a very thin 
wing at infinite Reynolds number, whereas the experiments 
were made with a 10-percent-thick wing at a Reynolds num- 
ber of approximately 5 X 10°. The effect of thickness and 
Reynolds number have not, as yet, been delineated for the 
oscillating case. 

As would be expected, it is shown theoretically that the 
larger the tunnel the less the effect of the walls. The critical 
frequency, however, is also reduced as the tunnel height is in- 
creased, but it is shown that tests may be made above the 
resonant frequency with no larger tunnel-wall effect than 
is found below the resonance. In addition, the range of 
influence of the resonant region is greatly reduced so that 
only a small range of frequencies need be avoided. Wall 
effects are shown theoretically to be more pronounced as the 
Mach number is increased and at high Mach numbers are 
found to be large even at frequencies well removed from 
resonance. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Ya., January 12 , 1055. 


APPENDIX 


REDUCTION OF INTEGRAL EQUATION TO THE CASE OF ZERO FREQUENCY 


Iii this appendix, the integral equation for the downwash 
for a wing oscillating in a compressible medium in the pres- 
ence of wind-tunnel walls is reduced to the zero-frequency 
condition. 

If equation (1) of the text is written as 

w(x)=lim -^2 f L(xo)[uK(M,z)-\ruK(M,z,ri)\dx 0 (Al) 
<o-*o pU J -i 

and the limit taken as co-^0, it will be found that all the 
terms of a>K(M,z) and uK(M,z,II) vanish except terms involv- 
ing Hi (2) • These terms become infinite; however, as co-^0, 
the asymptotic expansion for very small values of the 
argument may be used. Therefore, 
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The vertical induced velocity may then lie written as 
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Equation (A3) may be written as 


(A2) 


(A3) 


W ( x )=Y p r?H J , L ^ [ csch “ (l (5ll dx ° ( A4 ) 

Tlie additional induced velocity due to the presence of 
tunnel walls for the steady-state case in compressible flow 
is given by equation (40) of reference 13. Equation (A2) 


can be reduced to the same form by making the approxima- 
tion that the airfoil chord is small compared with the tunnel 
height. 
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